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We study a spin dynamics of two magnetic Chromium atoms trapped in a single site of a deep 
optical lattice in a resonant magnetic field. Dipole-dipole interactions couple spin degrees of freedom 
of the two particles to their motion in the site. The motion is quantized, therefore a trap geometry 
combined with two-body contact s-wave interactions influence a spin dynamics through the energy 
spectrum of the two atom system. Anharmonicity and anisotropy of the site results in a 'fine' 
structure of two body eigenenergies. The structure can be easily resolved by a weak magnetic 
dipole-dipole interactions. As an example we examine the effect of anharmonicity and anisotropy of 
the binding potential on the Einstein-de Haas effect. We show that the weak dipolar interactions 
provide a perfect tool for a precision spectroscopy of the energy spectrum of the interacting few 
particle system. 
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I. INTRODUCTION 



Dipolar magnetic interactions couple a spin of two par- 
ticles to their orbital motion in the trap [l| . A particular 
example of the spin dynamics driven by dipolar inter- 
actions and coupled to the orbital motion of the entire 
ferromagnetic sample is known as the Einstein-de Haas 
effect (EdH). The effect is a macroscopic illustration of 
the fact that spin contributes to the total angular mo- 
mentum of a system on the same footing as the orbital 
angular momentum. 

In the original experiment [H-Q a rotation of a ferro- 
magnetic rode was observed when the relative orientation 
of a magnetic field (forcing a magnetization of the sample 
along the field axis) and a magnetic moment of the sam- 
ple had been inverted. The observed rotation results form 
the fact that atomic spins adjust to a new orientation of 
the magnetic field. Conservation of the z-componcnt of 
the total angular momentum requires the orbital angular 
momentum of the sample to compensate for the change 
of the spin orientation. Therefore the orbital angular 
momentum has to depart from its initial zero value. The 
EdH effect strongly relays on the dipole-dipole interac- 
tions. These interactions, due to their anisotropy do not 
have to conserve the orbital angular momentum. Only in 
the presence of dipolar interactions initially non-rotating 
system can start to rotate. 

The EdH effect was also predicted for ultra cold 
trapped atomic gases [l-Q. Here it serves as a spectacu- 
lar manifestation of dipole-dipole interactions of atomic 
magnetic moments. This leads to a nontrivial dynam- 
ics where orbital motion is coupled to the atomic spin 
what allows to study many features of quantum mag- 
netism. Spin dynamics due to dipole-dipole interaction 
is inevitably associated with excitation of orbital mo- 
tion of an atom in the trap. The motion is quantized. 



there is therefore an energy gap separating initial and fi- 
nal states of different magnetization. The quantization 
of orbital energy of atoms significantly effects a spin dy- 
namics. This has been noticed in Q. 

Two scenarios are possible depending on the ratio of 
the energy gap to the dipolar energy. For a relatively 
dense gas of polarized Chromium atoms (10^^cm~^) ini- 
tially in S = 3,7715 ~ 3 state and for relatively weak 
traps (oj ~ 27r X lOOHz), the dipolar energy can exceed 
the excitation gap what makes the spin dynamics possi- 
ble. The atoms arc transfcrcd first to ms = 2 state and 
then, due to higher order processes and contact interac- 
tions, to other magnetic components [5[. 

Situation is different in a case of optical traps where ex- 
citation energies can be much larger then those in optical 
traps, or for species with small magnetic moments. The 
example is a ferromagnetic Rubidium whose magnetic 
moment is 12 times smaller then those of Chromium. In 
both above mentioned situations the dipolar energy is 
much less then the orbital excitation energy and the spin 
dynamics is suppressed due to the energy conservation. 

This obstacle can be overcame. As shown in the 
energies of both involved states can be tuned due to a 
Zeeman effect. The atoms prepared in = 3 state can 
be transferred to initially unoccupied ms = 2 compo- 
nents on the expense of the Zeeman, instead of the dipo- 
lar energy. In axially symmetric traps the phase of the 
ms = 2 atoms' wavefunction displays a characteristics 
vortex-like defect signifying orbital angular momentum 
in a final state. This scenario has been studied in the 
mean field limit, for a three component Rubidium spinor 
condensate Q. It has been shown that different final 
states, all involving vortex-like excitation, can be popu- 
lated via dipolar processes of the EdH type. Even highly 
excited trap states can be reached if energy and angu- 
lar momentum conservation are met by an appropriate 
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choice of the resonant magnetic field. 

Energy conservation makes the weak dipolar interac- 
tions highly sensitive, selective and tunable by means of 
external magnetic field. This requires however, a very 
high precision in control of both the value and the di- 
rection of the field applied. Stabilization of ultra weak 
magnetic field with precision below 1 mG is required for 
Chromium atoms [9[. The proposal exists suggesting 
that time dependent magnetic fields might relax these 
constrains [10| . Resonant charcter of demagnetization of 
initially polarized Chromim atoms in an optical lattice 
has been demonstrated recently [l^l • 

Recently, the EdH effect was suggested as a mecha- 
nism to create excited orbital states in the optical lattice 
PH . The mechanism leads also to an interesting physics 
described by Bose-Hubbard models with two spin species 
occupying different Wannier states. The axial symmetry 
of a single lattice site was assumed. This assumption 
seems to be quite obvious and is frequently used even in 
more 'exotic' situations [l3l |. 

In this paper we want to reexamine carefully the role 
of harmonic approximation and show its consequences 
on the spin dynamics, in particular for situations when 
higher orbital Wannier states are involved. We study 
elementary processes driven by dipolar interactions with 
two atoms trapped in a single site of the high square 
optical lattice. Such a scenario was considered in 
where it has been shown that 100% efficient resonant 
transfer of atoms to the orbital state is possible. We 
show that anharmonicity of the binding potential plays 
a very important role. In the context of ultracold atoms 
the effect of anharmonicity was discussed previously in 
p^ . Even small anharmonicity significantly affects the 
spin dynamics in a single site of an optical lattice. 

The paper is organized as follows. In section HIl we de- 
fine the Hamiltonian of the system as well as describe 
all the approximations and model parameters. Then in 
section Hill we discuss a resonant dipolar coupling allow- 
ing for efficient spin transfer in the channel where spin 
projection of one atom changes by one quantum. In sec- 
tion IIVI we present a detailed analysis of processes when 
spins of both interacting atoms change by two quanta. 
In particular we take into account separately the effect 
of anharmonicity of the trapping potential and, finally, 
both effect of anisotropy and anharmonicity of a single 
site of a square optical lattice. Final remarks are pre- 
sented in section IVl 



II. THE MODEL 

A. The Hamiltonian 

We study a simple system of two Chromium ^^Cr 
atoms in the S = 3 ground state in a single site of a 
2D square optical lattice created by a two pair of counter 
propagating laser beams of wavelength A and momentum 
k = 27r/A. A confinement in the z-dircction is provided 



by a harmonic potential of frequency fl^. We assume 
that the lattice depth is so big that we can neglect tun- 
neling to the neighboring sites. Through the paper we 
use the lattice units, i.e. the distance is measured in 
1/k = X/2tt, all energies in units of the recoil energy 
Eji = (ftfc)^/(2M), where M is the atomic mass. The 
units of time and frequency are h/ and Eji/h rcsspec- 
tively. The magnetic field unit is Bq ~ En/lgfis), where 
/^B is the Bohr magneton and g = 2 is a Lande factor. 
The Hamiltonian of the system has the following form: 

2 

H = ^ [-V^ + Vb(sin2 X, + sin^ y,) + n^zf)] 

i=l 

2S S 

+ 6{r,-r2)Y,gs ^ \S,M){S,M\ 

5=0 M = -S 

-{Si + S2)-B 

\ri - 

where = {xi,yi,Zi) and S{i) arc positions and spin 
operators of the i-th atom respectively. By n = r /r we 
denote the unit vector in the direction of r = ri — r2 = 
(x, y, z). Spin operators are represented by the spin 7x7 
matrices. To describe the contact two-body interactions, 
which depend on the total spin of both atoms, we use the 
molecular basis in the spin space: 

\S,M)^ (S,7m;S,m2\S,M)\S,m,)\S,m2), (2) 

mi ,7712 

where \S,Ai) is a two particle spin state with 
the total spin S and projection Ai. The symbol 
{S^rrii; S,m2\S , M) is the Clebsch-Gordan coefficient 
and \S, rrii) describes the state of the single atom of spin 
S and its projection to;. Note that due to the symmetry 
the total spin of the doublon, S can take only even values, 
S = 0,2,4,6. The parameter k is: k = {hn^)/{2ER). 
The strength of the contact interactions gs is propor- 
tional to the s-wave scattering length as in the given 
5-channcl: 

gs = STrkas- (3) 

Finally the dimensionless dipole-dipole interaction 
strength reads cf = fio/{4:TT){gfiB)^{k^/ER), where 
is the magnetic constant. 

The physical meaning of all the terms in the Hamil- 
tonian Eq.((T]) is straightforward. The first line corre- 
sponds to the single particle energy, Hq, including both 
kinetic and potential energies, while the second line de- 
scribes the contact interactions between particles. He- 
In principle, the contact interactions can lead to a spin 
dynamics. In this paper we assume that initially both 
atoms are prepared in the ground state of the lattice site 
in the presence of sufficiently strong magnetic field to 
assure a ferromagnetic polarization of the system. The 
initial spin state of the doublon is \S,Ai) — |6, 6). In 
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such a situation the spin dynamics is driven, at least 
initiahy, by the dipolar interactions only. The external 
magnetic field, the third line of Eq.([T]), describes the lin- 
ear Zeeman shift of single particles energies, Hz- The 
magnetic field allows for tuning of relative energies of the 
two atom system. In this paper we assume that magnetic 
field is align along vertical z-direction, perpendicular to 
the 2D lattice plane. Finally, the last line of Eq.([T]) cor- 
responds to the energy of the dipole-dipole interactions 
of both atoms, Hjj. The total Hamiltonian is the sum 
H ~ Hq + He + Hz + Ho ■ 

We examine the dipolar terms in more details stress- 
ing their role on the spin dynamics. The elementary pro- 
cesses driven by Hamiltonian Hd can be casted into three 
classes. 

1. Collisions not changing orbital angular momentum 
of the doublon, ACz = 0, where is the z-component 
of the total orbital angular momentum of the pair. These 
collisions are governed by the term whose spatial depen- 
dence has the form: 



ha 



1 - 



32^ 



(4) 



where and r = \r\. Note that this term has to be mul- 
tiplied by the appropriate operator acting in the spinor 
space: 

so = 5,(1)5,(2) - i (5+(l)5_(2) + 5_(l)5+(2)) . (5) 

The S±{i) = Sx{i) ± iSy{i) operator rises (-I-) or low- 
ers (— ) the magnetic quantum number of the i-ih 
atom by one quantum. The rising and lowering oper- 
ators fulfill the commutation relations for the spin, i.e 
[5+(fc), 5_(j)] = 2iSz{j)Skj- Note, that the spinor term 
So does not change the magnetization of the doublon, 
AA^ = 0. The dipolar interaction with ACz = are: 



H 



do 



ho ■ So 



(6) 



H^o evidently conserves the total angular momentum 
along the z-axis, A{£z + -M) = 0. Using explicit form 
of the spin matrices one can find that the strength of 
this interactions for both atoms in m = 3 components is 
- S^d^ = 3^d\ 

2. Collisions changing the z-componcnt of the or- 
bital angular momentum of the doublon by one quantum, 
ACz = ±1. These processes originate from the term of 
the following spatial dependence: 



h±i = - 



3d^ z{x ± iy) 



(7) 



The upper sign corresponds to ACz ~ 1, while the lower 
sign to ACz = ^1- This spatial term, h±i, has to be 
assisted by the corresponding spinor part 



Sti =•5.(1)5^(2) + 5^(1)5,(2). 



(8) 



The upper sign corresponds to processes with AA^ = 
— 1 , while the lower sign to the processes with AA4 = 



+1. The dipolar Hamiltonian describing processes with 
ACz = ±1 is a product of the spatial and the spinor 
terms: 



Hdi = h+i ■ .s_i + • s+i. 



(9) 



and clearly A(£, + M) ~ 0. The strength of these in- 
teractions for dipolar processes changing one atomic spin 
from TO = 3 and to = 2 is ~ = 3^/^d^. 

3. Collisions changing the z-component of the or- 
bital angular momentum of the doublon by two quanta, 
ACz = ±2 . These process originate from the term of 
the following spatial dependence: 



h±2 



3d^ {x ± iyf 



(10) 



which is accompanied by the spin term which changes 
the spin projection of the doublon by AA^ = ^2 



St2 



5^(1)5^(2). 



(11) 



The dipolar Hamiltonian describing processes with 
ACz ~ ±2 is a product of the spatial and the spinor 
terms 



H. 



d2 



h+2 ■ S-2 + h^2 ■ S 



+2- 



(12) 



Clearly the z-component of total angular momentum of 
the doublon A{Cz + M) = is conserved. The strength 
of this interactions for dipolar processes changing spins 
of both atoms from to = 3 to to = 2 is ^ Sd^ = 3d^ . 

This way we decomposed the dipolar interactions into 
different terms: 



Hd — Hdo + Hdi + i?, 



d2, 



(13) 



depending on the effect each term has on the orbital (or 
equivalently the spin) dynamics of the atom pair. 



B. General considerations 

To start with we assume that the external potential 
near the center of the site at r = can be approximated 
by the harmonic one. We keep only the lowest order 
terms in the expansion of the potential in Eq.([T]): 



Vtr = Vo(x' + y') + K^?? 



(14) 



The potential is both harmonic and axially symmetric 
with effective frequencies ~ '^x ^ ~ '^VVo and 



Viz = 2k. The aspect ratio is = k/vV^. 

In the following parts of the paper we assume that 
initially both atoms are prepared in the ground state of 
the external trap in the two body \S = 6,M = 6) spin 
state, what corresponds to a fully polarized system with 
both atoms in m = 3 state. The energy of the doublon is 
equall to sum of single partilce energies and contact inter- 
action. The initial orbital state 'I'ini(ri, r2) = (ri,r2|ini) 
docs not depend on angles of the vectors ri and r2: 

^'ini(ri,r2) = (l>Qo{xi,yi)(po{zi)(t)oQ{x2,y2)'Po{z2), (15) 
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where 0oo(a^, y)(po{z) is the ground state of the harmonic 
trap. In order to stress that axial wavefunction corre- 
sponds to a different trap frequency we separated out 
the axial term. 

To analyze further approximations we use the realis- 
tic parameters related to experiments of @. We assume 
that the optical lattice is created by the laser light of the 
wavelength of 532nm, which gives the unit of distance to 
be 84,67nm. For chromium ^^Cr {M ~ 87 x 10"27]jg) 
the recoil energy is En/h ~ 27r x 13, 5kHz. Typical 
height of the potential barrier is Vb = 25£'h what corre- 
sponds to Vo/fi = Stt X 336kHz for which the character- 
istic trap frequency Q.±_ = lOEn/h = 2tt x 135kHz. We 
assume in further discussion that the axial frequency is 
riz = IGEr/H =2nx 218KHZ. 

These parameters lead to the characteristic size of the 
system of the order of the harmonic oscillator length 
a± = 38nm and = 30nm. With two atoms per 
site the atom density is about uq ~ 2,8 x lO^^cm""^. 
Corresponding contact interaction in = 6 channel 
is Uc/^ ~ 27r X 18.9kHz while the dipolar interaction 
Ed /h ~ (fS - 27r X 0.37kHz. 

The largest energy is the separation between single par- 
ticles levels in the trap, i.e. Etr/h = ^± = 27r x 135kHz, 
The contact interaction energy is C/c/^ = 27r x 18.9kHz, 
and dipolar energy E^/fi ^ 2tt x 0.37kHz. The dipolar 
energy is the smallest energy scale. Thus we can ignore 
the term of dipolar interactions which does not change a 
total magnetization of both colliding atoms, i.e. we ne- 
glect processes with AAiz = 0. The term leads only to 
small, negligible corrections to contact interactions. 

Clearly, dipolar interactions leading to a spin dynam- 
ics assisted by an inevitable transfer of atoms to excited 
states of the trap are suppressed by the energy conser- 
vation because hCl± » Eo- The spin dynamics acti- 
vates if the energy gap is compensated by the Zecman 
energy. To shift relative energies of to = 3 and to = 2 
components by the amount of Ml±^ the magnetic field 
equal to Btr = 48.2mG is needed. Similarly, the Zee- 
man energy compensating for the contact interactions is 
Be ~ 6.75mG, while the dipolar energy corresponds to 
Bd ^ 25.7/iG. An extreme precise control of both di- 
rection and value of the magnetic field is necessary to 
control the spin dynamics in a trap. This makes a great 
challenge for experiments. 

For the initial axially symmetric ground state $ini = 
(ri,r2|ini) the processes with AAiz = —1 and AA4z = 
—2 lead to two different orbital states of the doublon. In 
the lowest order of the perturbation in dipolar interac- 
tions, the angular dependence of the two atoms' relative 
wavefunction (coupled to the ground state by dipolar in- 
teractions) is determined by the angular dependence of 
the dipolar terms, Eqs.Q, (fTOj) . The state excited due 
to AMz = — 1 dipolar process has energy of the order 
of + flz while the energy of the state reached by the 
system via the second channel, AAiz = —2, is 2^1 ±. If 
(flz — Sl±) is much larger than dipolar energy then by 
matching the resonance condition one can activate the 



first or the second channel separately (not both of them) . 
Therefore we will analyze the process with AAi^ = — 1 
and AAiz = — 2 separately. 

Our general strategy is to find the lowest energy state 
which is coupled by dipole-dipole interactions to the |ini) 
state as it follows directly from the structure of the dipo- 
lar Hamiltonian. At this step we determine a single parti- 
cle basis which is needed to describe all relevant two-body 
states. Then we diagonalize the full Hamiltonian H in 
the chosen basis and find cigenenergies and eigenvectors 
of the doublon. At every resonance there exists an eigen- 
vector of the Hamiltonian which is an equal-amplitude 
superposition of the initial and the final state. The two 
states forming the superposition are coupled by the dipo- 
lar interaction. 

Alternatively, the physical insight can be gained with 
the help of a perturbative analysis. The state determined 
by the structure of dipole-dipole term is not necessar- 
ily the eigenstate of the single particle and the contact 
Hamiltonian. The goal is to decompose this state in the 
eigenstate basis of Hq + Hz + He- A number of dipo- 
lar resonances equals to a number of different cigcnstates 
involved in the decomposition. Each component has dif- 
ferent eigenenergy and therefore becomes resonant at dif- 
ferent magnetic field. 

To the 'zero order' approximation the energy of a final 
state - the state accessible by the dipolar interactions - is 
given by the sum of the single particle energies (f2_L -I- ^Iz 
in the AM.z = — 1 channel or 20^ in AMz = —2 chan- 
nel) . The single particle energy gives an estimation of the 
resonant magnetic field. However, contact interactions 
adds important contributions to the energy of the two 
atom system. This effect results in splitting of energies 
of different components of the final state and multiple 
resonances do appear. 

The strength of the resonant dipole-dipole interactions 
can be characterized by the self energy of the initial |ini) 
state of energy Eq{B): 

Esf{B)^J2 , ^' . (16) 

l^'i) and Ei{B) are 'bare' cigcnstates and cigenener- 
gies of the two-body Hamiltonian without dipolar in- 
teractions, {Hn + He + -ffz)|*i) = E,{B)\^i), while 
Di = (ini|ij£)|^'i) is the dipole-dipole interaction ma- 
trix element. The self energy contains contribution from 
a virtual transfer of the doublon from |ini) to the inter- 
mediate l^'i) state. Evidently, this process contributes 
only at the vicinity of the resonance, where: 

E,{Bres) - EoiBres) = D,, (17) 

The particular form of the self energy results from the 
fact, that at every resonance we deal with a two level 
system, and then the denominator in Ea. p6p gives the 
energy shift of the |ini) state. 



III. DIPOLAR PROCESSES WITH AMz = -1 

The spatial part of the state of the doubloii coupled 
by the dipole-dipole Hamiltonian Eq.([7]) in AM^ = ^ 1 
channel to the initial state |ini) is z(x + iy), where 
X, y, z, are relative coordinates of two atoms. This state 
is created by a two-particle operator: 



(18) 



acting on the particle vacuum. The bosonic operators 
in the above formula create a particle in one of the fol- 
lowing harmonic oscillator states: creates a particle 
in the ground state (j)oo{x,y)(po{z), pi and pj^ create an 
atom in one of the two p-shell states - (j)ioix, y)ipo{z) and 
(j>oi{x,y)Lpo{z), while and d^^ create a particle in the 
d-shell states - (l>io{x,y)(pi{z) or (j)oi{x,y)(pi{z). ipi{z) is 
the first excited state in the axial direction. The opera- 
tor vl allows for a compact representation of the 'vortex' 
state with Cz = 1. The vortex is created in a relative co- 
ordinate space of the doublon. The first term in Eq. pB]) 
describes the state in which one atom takes one excitation 
quanta in the axial direction, while the second carries the 
vortex like excitation. The second term corresponds to 
a process in which one atom experiences both types of 
excitations while the second is a spectator (at least as 
the spatial wavefunction is considered). 

The operator v\ defines only a spatial part of £, = 1 
vortex state: 



\Ci)=vl\0) 



(19) 



To account for a total state of the system including spin 
degrees of freedom, the |£i) state has to be multiplied by 
a totally symmetric spinor component with one atom in 
\m = 3) and the other in \m = 2) spin states. The final 
state of the system is: 



1 
72 



(|3)|2) + |2)|3))|A 



(20) 



The state \Ci) is a sum of two components. Evidently 
in the harmonic trap the single particle energies of both 
components are the same. The most important observa- 
tion is that contact interaction vanishes in the |£i) state 
because its wavefunction vanishes identically if ri = r2, 
i.e. (r, r|£i) = 0. This state is not coupled by con- 
tact interaction to any other state of the doublon in the 
trap. It is therefore the eigenstate of both the single 
particle and the contact Hamiltonians. Dipole-dipole in- 
teractions couple efficiently the state to the initial one if 
energies of the involved states do match. This can hap- 
pen if = ^± + + {2 + 3)B equals to the energy of 
the initial state Eini = U33 + 2{3B). The resonance con- 
dition is met at magnetic field Bres = i^z + ^^_l) — £^33- 
t/33 ~ {gQ/8)yy^l\ftz/Tr^ is the contact interaction of two 
particles in m = 3 state and the ground state of the har- 
monic potential. 




FIG. 1: Dipolar resonances with AA4z ~ —1 in the harmonic 
trap. The excitation energy of the two particle system in 
the harmonic trap as a function of external magnetic field. 
The horizontal line corresponds to the |ini) state with both 
atoms fully polarized |3)|3). At zero magnetic field the state 
becomes the ground state for a given magnetization. Only 
the lowest energy state \£i}, the one which is coupled to the 
|ini) state by AAlz = — 1 dipole-dipole interaction term is 
shown (left upper panel). Notice avoided crossing as well 
as a maximum of the self-energy curve (right upper panel) 
signifying strong coupling of both involved states at resonant 
magnetic field. Lower panel: reduced single particle density 
matrix component (the diagonal element with magnetization 
m = 2) of the state |£i); left panel - density cut at xy plane 
for 2 = right panel - density at xz plane (j/ = 0). 



The state is also robust if anharmonicity of the lat- 
tice site potential is taken into account. In our case the 
lattice potential depends only on x and y coordinates. 
The z-direction separates. This fact has very important 
consequences: d^z and dyz states are separable. Their 
single particle energies are sum of p-band energies. For 
this reason the two components of the |£i) state, i.e. 
the one with two excited particle and the second with 
only one excited particle, are degenerate. Anharmonic- 
ity of the lattice potential results in an equal shifts of 
both components of the |£i) state and does not lift their 
degeneracy. Moreover, this conclusion holds even when 
axial symmetry of the external potential is substituted 
by the Z2 symmetry of the single site of the lattice. Z2 
symmetry is sufficient to assure degeneracy of Px and Py 
orbitals. 

The latter will not be true if Px and Py orbitals have 
different energies. In such a case the |£i) state will be 
a superposition of two eigenstates of different energies, 
and each component could be separately addressed in a 
dipolar process by an appropriate tuning of the magnetic 
field. Two resonances will be observed in such a case. 
The only effect of the lattice potential is to shift the value 
of the resonant magnetic field. 

Diagonalization of the full Hamiltonian in the oscilla- 
tory or Wannicr basis supports the above discussion. We 
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find only one low energy state coupled by dipole-dipole 
interactions to the initial one. In Fig [1] we show eigcnen- 
ergies of the two states. The avoided crossing, signifying 
the dipolar coupling is clearly visible. In the lower panel 
the single particle density in m = 2 spin component is 
shown {xy and xz cuts). Note the node in z = plane 
what is typical for the first excited state in axial direc- 
tion. The single particle density does no resemble the 
structure typical for a single particle vortex state. This 
is because the vortex is created in the relative coordinates 
of two particles. At the top on the right panel of FiglT] 
the self energy of the initial state is plotted. The hight of 
the peak visible at resonance position indicates the value 
of the dipolar matrix element coupling the states. If the 
magnetic field is tuned to the resonant value the system 
will perform Rabi oscillations between the initial and the 
final state. The final state reached in the dipole-dipole 
interaction is highly entangled, Eg. ([^0]) . 

IV. DIPOLAR PROCESSES WITH AM^ = -2 

A. Harmonic trap 

In this section we again assume that the doublon is 
trapped by the harmonic potential as given by Eq. lfT^ 
and we examine collisions which change the magnetiza- 
tion of the system by two quanta, AAi^ = —2. The 
two interacting particles go simultaneously from |3)|3) to 
|2) |2) spin state. Spatial part of the relative wavefunction 
in the lowest accessible final state must be proportional 
to ~ (x -|- iy)^ in order to compensate for a decrease 
of the spin of the doublon by two quanta. This is the 
wavefunction of a doubly charged vortex in the relative 
coordinates space. This simple state (normalized) if writ- 
ten in the harmonic basis looks quite complicated, so it is 
convenient to describe it using the second quantization 
picture. The 'vortex' state is created by the two-body 
operator u|: 




In addition to defined previously bosonic operators creat- 
ing particle in s and p shells we define d shell operators: 
dj,^, d^yy, and d^y creating a particle in the second ex- 
cited state of a harmonic oscillator, i.e. (j)2oix,y)ipo{z), 
4>02{x,y)ip^{0), and 0ii(a;, y).^o(^) respectively. 

The first line describes the state with one particle in 
the ground s state and the second particle in the excited 
d-shell state with Cz = 2. The second line corresponds 
to a state with two particles occupying the Px + ipy vor- 
tex. Let us note that introduced here states, abbreviated 
by s, 

Px, Py, dxx, dyy, and dxy, span a basis of orbital 
states which is sufficient to describe the lowest dipolar 
excitations in AAiz = — 2 channel. 



The orbital basis has to be completed by the spin 
states. These are the initial |3)|3) and the final |2)|2) 
states. The further spin dynamics due to the spin de- 
pendent contact interaction can subsequently induce a 
transition from |2)|2) to (|1)|3) + |3)|l))/v^ state. How- 
ever, the process can be suppressed due to the energy 
mismatch if to = 2 state is shifted by light. Therefore for 
a moment we neglect contact spin-changing collisions as 
it can be a realistic assumption. 

We diagonalize the total Hamiltonian Eq. ([T]) in the ba- 
sis of the two particle states spanned by the introduced 
above single particle s, p and d orbitals in the initial 
TO = 3 and the final to = 2 spin component for different 
values of external magnetic field. 

However to gain an insight into physical processes in- 
volved we shall describe the results using perturbative ap- 
proach with a dipole-dipole interaction energy as a small 
parameter. First we diagonalize the single particle and 
contact Hamiltonian Hq + Hz + He- As expected, the 
state |£2), defined as: 

l^2)=4|0), (22) 

is the eigenstate of Hq+Hz + Hc- Moreover, the contact 
energy vanishes in this state. This is a unique feature of 
the harmonic potential and has some important conse- 
quences on the spin dynamics as it will be shown later. 
The energy of this state — 2il^ + 2{2B) equals to 
the energy of the initial state Eini = C/33 + 2(3-6) at mag- 
netic field Bres = n^- C/33/2. C/33 = (56/8)v/^^IfV^ 
is the contact interaction of two particles in to = 3 state 
and the ground state of the harmonic potential. 

Evidently, by the construction, [£2) is coupled to the 
initial state: D = (£2|^+2|iiii) 7^ 0. By diagonalization 
of the full Hamiltonian with the dipole-dipole term in- 
cluded, we indeed find that [£2) state is the only one 
coupled to jini) state by the dipole-dipole interactions. 
The two coupled states, |ini) and |£2) states form an ef- 
fective two level system. The energies of the two states 
for different values of external magnetic field are shown 
in Figl2l At the resonant magnetic field, Bres, there is 
an avoided crossing. Of course, there are higher excited 
states, not included in our basis, which are coupled by 
dipolar interactions to the initial state. This coupling 
takes place, however at higher magnetic fields which we 
do not consider here. 

At the top right site of FigI5] we show the self energy 
of the |ini) as a function of the magnetic field. The self 
energy reaches a maximum at the resonance. The mag- 
nitude of Esf at the resonance is equal to D. Similarly, 
a width of the resonance is also equal to D. 

Weakness of the dipolar coupling allows for a precise 
tuning of the two atom system. Initial population of 
the TO = 3 state will oscillate periodically between two 
coupled states with a Rabi frequency proportional to the 
dipolar energy D. At half of the period all the population 
will be transferred to the to = 2 state, |£2)- The single 
particle density corresponding to this state is shown in 
the lowest panel of Fig 1^1 Note a shallow minimum at 
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FIG. 2: Dipolar resonances with AyVfz — —2 in the har- 
monic trap: excitation energy of the two particle system in 
the harmonic trap as a function of external magnetic field. 
The horizontal line corresponds to the |ini) state with both 
atoms fully polarized 13)|3). At zero magnetic field the state 
becomes the ground state for a given magnetization. Only 
the lowest energy \jC2) state coupled to the |ini) state by 
AA4z = — 2 dipole-dipole interaction term is shown (left up- 
per panel). Notice avoided crossing as well as a maximum 
of the self-energy curve (right upper panel) signifying strong 
coupling of both involved states at resonant magnetic field. 
Lower panel: reduced single particle density (in xy plane) of 
I £2) state coupled to the |ini) state. 



the trap center. The vortex \C2) state is a vortex in the 
relative coordinate space. This is no vortex in a position 
space of a single particle. The doubly occupied vortex 
state Px + ipy involves excitation of the center of mass of 
the doublon. This is however not possible with any two 
body interactions in the harmonic potential. 



Finally we return to the observation that the con- 
tact interaction in |£2)|2)|2) state vanishes, i.e. U22 = 
322('C2|i5(ri - r2)|£2) = 0, where 322 = (655 + ^g^)/!!- 
This is because the two body wavefunction of this state 
vanishes if ri = r2. The contact interaction leading to 
the transfer of one atom to m = 3 and the second one to 
m = 1 state, C/22,31 = 522,31 ('C2|<5(ri - r2)|£2) = 0, also 
vanishes, where 522,31 = ■\/30/ll(36 — 94)- Therefore, in 
the harmonic axially symmetric trap the spin dynamics 
of the |£2)|2)|2) state due to the contact spin changing 
collisions is forbidden. There is no need to detune m = 2 
state with respect to m = 3 and m = 1 states to prevent 
contact spin changing collisions. This fact is in a con- 
tradiction with experiments with Chromium atoms in an 
optical lattice and clearly indicates that harmonic ap- 
proximation does not describe correctly a spin dynamics 
in optical lattices. 



B. Anharmonic trap 

If in the expansion of the lattice site potential we keep 
the higher order terms then we get a potential which is 
neither harmonic nor axially symmetric: 

Vir) = Vop^ - ^Vop" + K^z^ + '^Vox'^y^ (23) 



We introduced notation: p = 



y 



The second term 



^ p"^ introduces anharmonicity while the last term brakes 
the axial symmetry. To separate effects of the anhar- 
monicity from effects of the anisotropy we neglect the 
last term and consider the model anharmonic potential 



lvop' + ^'z\ 



(24) 



The combined effect of anharmonicity and anisotropy of 
the square lattice site will be discussed in the subsequent 
section. 

A measure of anharmonicity is given by a ratio of the 
anharmonic to the harmonic energy at a typical distance 
PO: 7 = jPo- The radial frequency of the harmonic term 
fl± = 2\/Vq defines the spatial extension of the ground 
state wavefunction pq = 1/V^-L- Therefore the anhar- 
monicity coefficient is: 



1 1 



7 = 7 



6 



(25) 



Obviously anharmoniciy 7 is small for deep optical lat- 
tices. The energy levels of anharmonic potential are not 
equally spaced. For this reason the state |£2) is no longer 
the cigenstate of the Hamiltonian Hq + Hz + He ■ 

To construct two-body cigcnstates we find numerically 
all the s, p and d-shell single particle eigenstates of the 
anharmonic potential. Then the full two-body Hamilto- 
nian is diagonalized in this basis. There are two states 
|uai) or jwaa) which are coupled to the initial one at two 
different magnetic fields. The resonances are visible as 
avoided crossing in Fig (31 At each resonance the system 
is effectively a two level one - the population will oscil- 
late between the initial state, |ini), and the resonantly 
coupled one, i.e |uai) or Iwaj). 

The single particle densities of the two resonantly cou- 
pled states for Vb = 25£^fl are presented in the lowest 
panel of FiglS] Let us note, that one particle density of 
the \va2) state resembles a singly quantized vortex, but 
does not reach the zero value at the center. In the middle 
panel of Fig 13] we show the self energy of the |ini) as a 
function of the magnetic field. The self energy reaches a 
maximum at every resonance. 

The presented here numerical results can be under- 
stood with the help of the perturbative approach with a 
dipolar energy being a small parameter. The state I/I2) 
can be viewed as the superposition of the two states: 



l>C2) 



1 



V2 



(26) 
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FIG. 3: Dipolar resonances with AtMz = —2 in the an- 
harmonic trap: excitation energy of the two particle system 
in the anharmonic trap as a function of external magnetic 
field. The horizontal line corresponds to the |ini) state with 
both atoms fully polarized |3)|3). At zero magnetic field the 
state becomes the ground state for a given magnetization. 
Only the lowest energy states coupled to the |ini) state by 
AA^z — —2 dipole-dipole interaction term are shown (upper 
panel). Notice two avoided crossing as well as two maxima 
of the self-energy curve (middle panel) signifying strong cou- 
pling at resonant magnetic field. Lower panel: reduced single 
particle densities (in xy plane) of the final states accessible 
via AA4z = —2 resonant dipolar processes corresponding to 
the two resonant values of the magnetic fields as indicated in 
the middle panel. 



where d-shell state 



Vd) is: 



\vd) = i (fit (4, + - 4,) ) (27) 

and p-shell state \vp) has the following form: 

■ ■ ^ ((pt)2+,2ptpt_(pt)2)|o), (28) 



2\/2 



In formula Eq. ((27|) and Eq. p8|) operators s , p , d are not 
the same as those in the harmonic case however we use 
the same notation. Now the operators create a particle 
in s , p , d states of the anharmonic potencial. 

The states have different single particle energies - d 
and p states arc shifted differently due to anharmonicity. 
Both states arc the cigenstates of the z-component of the 
orbital angular momentum of two particles. \vd) is the 
state with one particle in the ground and the second in 
the excited d-state with >Cj; = 2. \vp) is the state with 
two particles occupying the same p orbital with £2 = 1. 
Evidently the single particle density of the second state is 



typical for the singly charged vortex state in the harmonic 
trap. 

In a harmonic trap the two states have the same en- 
ergies. Their relative shift introduced by anhamonicity, 
5E, can be estimated using the perturbative approach 
with 7 being the small parameter: 



SE = 27f7j 



(29) 



The energy of \vd) state is smaller then the energy of \vp) 
state. 

These states are dressed by the contact interaction 
whose strength is: 



^22 - ^§^v/^^^J 



(30) 



where 322 = (6g6 + 5g4)/ll. The dressed eigenvectors 
and eigenenergies of this effective two level system are 
solutions of a simple eigenvalue problem: 



-SE L/22 \ / aa 



U22 



Pa 



(31) 



The two eigenstates 

Pa\vd) 

states 



aaV'd) ~ Pa\vp) and l^aa) 



+ Cia\vp) are superpositions of the two coupled 



and 



If the contact interaction significantly exceeds the en- 
ergy splitting U22 >> SE/2, i.e. for: 



Vo » 



256 



9 {922rn/' 



(32) 



then harmonicity is recovered in a deep optical lattice. 
Note the crucial role of the axial trapping frequency. In 
a prolate geometry harmonicity is recovered for much 
deeper traps than in a case of an oblate geometry. For 
Clz = 16 as used here, the anharmonicity of the opti- 
cal lattices becomes negligible if Vb >> 29Er. In this 
limit the eigenvectors of the Hamiltonian are |uai) — 
i\vd) - \vp))/V2 and \va,) ^ {\vd) + \vp))/V2. Note, that 
the eigenstate jwai) corresponds to relative coordinates 
excitations and is the same as the state reached by the 
dipolar interactions in the harmonic trap, l^ai) — \^2)- 
The second one, l^aj), decouples from the initial |ini) 
state as corresponding dipolar matrix element vanishes. 
This states corresponds to the center of mass excita- 
tion. Obviouusly excitation of the center of mass is not 
possible in a harmonic trap by any two-body interac- 
tions. We checked that even for deep optical lattices, 
i.e. Vq ~ 4QEji, the effect of anharmonicity cannot be 
neglected. 

Contact interaction becomes negligible if, U22 << 
SE/2. Then the eigenstates of the interacting systems 
can be approximated by the jwa^) ~ \vd) and \va2) — \vp)- 
The one-particle density of the latter state is characteris- 
tic for a singly quantized vortex - it is therefore a perfect 
candidate to observe the EdH effect. The case shown in 
FiglSl corresponds to a crossover situation SE ~ 1/22- 

Till now we did not take into account the spin dynam- 
ics due to the contact interaction in the final |2)|2) spin 
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space of the doublon. This simphfication is justified if 
contact interaction couphng the two spin subspaces are 
'turned off' by a Ught shift of the rn = 2 state. If it is 
not the case then the contact interaction foUowing the 
dipolar process can further transfer the doublon to |3)|1) 
spin state. 

Oversimplified picture of |2)|2) to |3)|1) coupling is 
based on the observation that the eigcnstates |uai) and 
jwaj) of two body contact Hamiltonian have their anal- 
ogons of the same eigenenergy and the same wavefunc- 
tions in |1)|3) subspace. If so, then the state |wai)|2)|2) 
{i = 1,2) is coupled to |uaj(|l>|3) + |3)|1))/V2 state due 
to the contact term 522,31 (woi |<5(ri — r2)||t'ai). Then the 
form of the two-body wavefunction in |3)|1) space is in- 
herited from the |2)|2) subspace. 

Situation is more complicated however. The con- 
tact interactions are different in different spin subspaces. 
They are proportional to 1722 = (6(76 + 5(74)/ll or to 
531 = (556 + 654)711 in |2)|2) or in |3)|1) spin sub- 
space respectively. For this reason the vectors diago- 
nalizing contact interactions in |3)|1) space are different 
then |wai) and [waa) diagonalizing the contact interac- 
tions in the |2)|2) space. Their eigenenergies are also 
different. Note however that difference in the contact in- 
teraction strength is 522—531 = 522,3i/'\/30, i.e. by factor 
y/SO smaller than the term coupling both spin subspaces 
522.31- For this reason, the simplified discussion is not 
totally wrong. 

Exact treatment requires taking into account the two 
pairs of vectors from different subspaces. In order to find 
the states of the doublon dressed by the contact inter- 
actions one has to diagonalize the 4x4 matrix in the 
basis of 1?;^) and \vp) states in the two involved spin sub- 
spaces. The splitting of energies of the four eigenstates 
is determined by the contact interactions. The eigen- 
states are superpositions of different spin states of the 
doublon. Each of them can be reached in the dipolar re- 
laxation process if the external magnetic field is tuned to 
the given resonant value. 

In FiglH the four resonances are clearly visible for 
Vq = 25Efi. The single particle densities correspond- 
ing to the every resonant state of the doublon are shown 
in the lowest panel. Let us observe first that parameters 
we choose are in the crossover regime. The anharmonic 
shift of energies, although does not dominate the con- 
tact interactions, is still important. The low magnetic 
field resonant state, in both spin subspaces, is dominated 
by the d-shell \vd) orbital with one particle in the trap 
ground state (low panel of FigU]). That is why the single 
particle density is large in the trap center. The higher 
field resonant state is dominated by the p-shell excita- 
tion. This is the reason of the density minimum at the 
center. 

A small anharmonicity of the binding potential signifi- 
cantly modifies the spin dynamics. In the harmonic trap 
the two atoms oscillate between |ini)|3)|3) and |£2)|2)|2) 
state at the resonant magnetic field B^es- The I/I2) state 
is characterized by the orbital angular momentum Cz = 2 
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FIG. 4: Dipolar resonances with AA4z = —2 in the anhar- 
monic trap: excitation energy of the two particle system in 
the anharmonic trap as a function of external magnetic field. 
The horizontal line corresponds to the |ini) state with both 
atoms fully polarized |3)|3). At zero magnetic field the state 
becomes the ground state for a given magnetization. Only the 
lowest energy states coupled to the |ini) state by AA4z = — 2 
dipole-dipole interaction term and spin dependent contact in- 
teractions are shown (upper panel). Contact interaction in 
|2)|2) copling to l/y2(|3)|l) -|- |1)|3) spin state is included. 
Notice four avoided crossings as well as four maxima of the 
self-energy curve (middle panel) signifying strong coupling at 
resonant magnetic field. Lower panel: reduced single parti- 
cle densities (in xy plane) of the final states accessible via 
AMz = —2 resonant dipolar processes and the spin depen- 
dent contact interactions corresponding to the four resonant 
values of the magnetic fields as indicated in the middle panel. 



with no center of mass excitations. For a realistic lattice 
depth this single resonance is split into four resonances. 
The spatial wavefunctions of the doublon are, in all four 
cases, eigenstates of the Cz component of the orbital an- 
gular momentum (with the center of mass excitation). 
The states are entangled. This entanglement can be 
use full, particularly after applying a Stern-Gcrlach mag- 
netic field separating spatially different entangled mag- 
netic components of the doublon. 



C. Single site of a square lattice 

Finally we are ready to discuss a realistic situation of 
two atoms in a single site of the square optical lattice 
where both: the effect of anharmonicity and anisotropy 
are present simultaneously. Again, in order to find final 
states of doublon which can be reached by dipolar inter- 
actions from the initial |ini)|3)|3) state we have to find 
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the eigenstates of the single particle and contact Hamil- 
tonian Hq + Hz + He- To this end we have to adjust the 
single particle basis to the lattice case. Instead of oscil- 
latory wavefunctions we use the Wannier states (except 
of the z-component part which is the ground state wave- 
function of the harmonic oscillator in the z-direction. To 
simplify a notation we omit the z-dependence of the basis 
wavefunction - but we keep it in our calculations) . If we 
limit our basis to the three lowest bands in the lattice 
then the best approximation to the vortex state is the 
state created by the operator v\ as defined in Ea. ((2T|) . 
The operators s\ p^, p]y, dj,^, d)yy, dly have slightly dif- 
ferent meaning now. They generate a particle in the 
Wannier states Wiix) rather than in the harmonic os- 
cillator states. For example creates particle in the 
^-00 = Wo (a;) Wo (y) state, 4, in the *2o = m{x)Wo{y) 
state, or dl.y in the = Wi(a;)Wi(y) state, etc. 

It is convenient to rearrange the terms in the expression 
for the vortex state \C2) according to their single particle 
energies in the lattice, i.e: 



7.4 7.6 7 



2 8.4 8.6 8.1 



\C2) - |vi) + |V2) + |V3), 



where: 



1 



|V2> 
|V3) 



2^/2 
1 



(4.-4,))|o), 



^ {pIpI-pIpD |o), 

{Suiy-Plpl)\0). 



(33) 

(34) 
(35) 
(36) 



The single particle contribution to the energies of these 
states are Ei = Es + Ed, E2 = ^Ep, and £'3 = E2 respec- 
tively, where i?-^^^ = (Wi|i?o|Wj;) is the mean energy of 
particle in the given Wannier state. The single particle 
energies of the last two states are equal because of the Z2 
symmetry of the problem. Anharmonicity is responsible 
for the inequality Ei < i?2- 

The goal is, to find eigenstates of the contact Hamil- 
tonian and decompose \C2) state in the basis of dou- 
blon eigenstates. Note, that the two-body wavefunction 
(ri,r2|v3) vanishes identically if ri = r2 so does the con- 
tact interaction in this state as well as all possible con- 
tact interactions coupling this state to the other states. 
The state IV3) is therefore the eigenstate of the two-body 
Hamiltonian Hq + Hz + He ■ 

On the contrary contact interaction couples the vec- 
tors |vi) and |v2). The states |vi) and |v2) are analogous 
to \vd) and \vp) states discussed in the previous section. 
Their two linear superpositions diagonalize the contact 
Hamiltonian. The two eigenvectors can be found numer- 
ically. 

In the absence of coupling to the |3)|1) subspace, there 
are three low energy states which can be accessed due 
to the dipole-dipole interactions. The three two-body 
eigenenergies (meassured from the energy of the initial 
state) are shown in FiglS) Avoided crossings signifying 
dipole coupling are clearly visible. In the lower panel 
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FIG. 5: Dipolar resonances with AA4z = —2 in the single 
site of a square lattice: excitation energy of the two particle 
system in a single lattice site as a function of external mag- 
netic field. The horizontal line corresponds to the |ini) state 
with both atoms fully polarized |3)|3). At zero magnetic field 
the state becomes the ground state for a given magnetization. 
Only the lowest energy states coupled to the |ini) state by 
AMz = —2 dipole-dipole interaction term are shown (upper 
panel) . Notice three avoided crossing as well as three maxima 
of the self-energy curve (middle panel) signifying strong cou- 
pling at resonant magnetic field. Lower panel: reduced single 
particle densities (in xy plane) of the final states accessible 
via AMz ~ —2 resonant dipolar processes corresponding to 
the three resonant values of the magnetic fields as indicated 
in the middle panel. 



the single particle densities of the states coupled to the 
initial one are shown, while in the middle panel the self 
energy of the initial state is presented. The self energy 
depicts position, strength and width of every resonance. 
Note, that for axially symmetric harmonic trap, all three 
resonances will coincide forming Cz = 2 vortex state in 
the relative coordinates space of the two atoms. 

If we account for the contact interactions coupling |vi) 
and |v2) states in |2)|2) spin subspace to analogous states 
in |3)|1) subspace the additional channels of the spin dy- 
namics do appear. The situation is analogous to the one 
discussed in the case of anharmonic trap. Again, the 
4x4 eigenvalue problem leads to four different eigen- 
vectors with entangled spin and spatial wavefunctions. 
Together with |v3) state there arc five low energy dipolar 
states resonantly coupled to the initial one in the single 
site of optical lattice. State |v3) is not shifted by contact 
interactions, 'lives' entirely in the |2)|2) spin subspace. 
The resonance corresponding to this state is visible at the 
top of FiglS] The other resonances correspond to dipolar 
coupling of the initial state |ini) to the four eigenstates 
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FIG. 6: Dipolar resonances with AA4z ~ —2 in a single 
site of a square lattice: excitation energy of the two par- 
ticle system in a single lattic site as a function of exter- 
nal magnetic field. Contact interactions in |2)|2) copling to 
l/\/2(|3)|l) -|- 11)|3) spin state is included. The horizontal line 
corresponds to the |ini) state with both atoms fully polarized 
1 3) 1 3). At zero magnetic field the state becomes the ground 
state for a given magnetization. Only the lowest energy states 
coupled to the |ini) state by AMz = — 2 dipole-dipole in- 
teraction term and spin dependent contact interactions are 
shown (upper panel). Contact interaction in |2)|2) copling 
to 1/^(|3)|1) -I- |1)|3) spin state is included. Notice five 
avoided crossings as well as five maxima of the self-energy 
curve (middle panel) signifying strong coupling at resonant 
magnetic field. Lower panel: reduced single particle densities 
(in xy plane) of the final states accessible via AMz ~ —2 
resonant dipolar processes and the spin dependent contact 
interactions corresponding to the five resonant values of the 
magnetic fields as indicated in the middle panel. 



of the Hamiltonian Hq + He + Hz, obtained by numer- 
ical diagonalization in the Wannier basis, in Figl6l The 
single particle densities of the states accessible by the 



dipole-dipole interactions are shown in the lower panel. 
The single particle densities depend on the z-componcnt 
of the spin of the atom. 



V. FINAL REMARKS 

We have shown that dipole-dipole interactions are 
very effective probe of the two-body states of Chromium 
atoms in the external trapping potential. They are very 
sensitive to anharmonicity and symmetry of the single 
site. The anharmonicity of the single trap becomes very 
important in particular when higher Wannier states with 
non vanishing angular momentum are concerned. The 
contact interaction between the degenerate states in the 
site lifts this degeneracy and introduces a fine structure 
of the energy spectrum of the system. Spectroscopy of 
this fine structure can be done with the help of dipo- 
lar resonances. Observation of spin dynamics driven by 
contact interactions is a precise method of determination 
of contact interaction for both Fermi and Bose atoms 
[l6l . [TtI . We show that spin dynamics due to dipolar 
interactions at resonant magnetic field is yet another 
perfect tool for precise determination of the scattering 
lengths of involved atoms. This can be very important 
for measurement of various scattering length for dipolar 
elements such as Dysprosium or Erbium [3, • Gener- 
alization of the approach presented here to more atoms 
per site or Fermi species is straightforward. 
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